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SUMMARY 

The  mathematical  description  of  the  geometry  for  forward  specular 
scattering  of  radio  waves  by  field-aligned  ionization  irregularities 
has  been  extended  to  cover  transequatorial  ray  paths.  The  equations 
presented  include  expressions  for  paths  involving  either  one 
scattering  point  in  the  same  geomagnetic  hemisphere  as  the  trans¬ 
mitter,  or  two  scattering  points,  one  on  each  side  of  the  geomagnetic 
equator . 
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1 .  INTRODUCTION 

Magnetic  field-aligned  ionization  has  been  detected  in  the  ionosphere  for  many 
years  and  by  several  experimental  techniques.  Such  irregularities  are  most 
prevalent  in  the  auroral  region  but  have  also  been  found  to  exist  at  low 
latitudes  (refs.  1,  2).  The  generally  accepted  form  of  the  ionization 
irregularities  is  of  a  tube  or  column  of  increased  or  depleted  electron 
density  confined  and  elongated  along  a  magnetic  field  line.  Radio  waves 
incident  on  magnetic  field-aligned  ionization  in  a  direction  perpendicular  to 
the  magnetic  field  lines  are  scattered  back  in  the  direction  of  incidence,  a 
property  which  has  formed  the  basis  for  most  experimental  investigation  of  the 
phenomenon . 

In  addition  to  this  aspect-sensitive  scattering,  incident  radiation  can  undergo 
forward  scattering  when  the  direction  of  propagation  makes  an  angle  of  other 
than  90°  with  the  magnetic  field  lines.  Experimental  observations  have 
demonstrated  that  forward  scattering  is  specular  or  nearly  specular,  i.e.  the 
angle  of  scatter  (or  reflection)  is  equal  to  the  angle  of  incidence.  Thus  the 
scattered  signal  describes  a  conical  surface  surrounding  the  column  of 
ionization,  which  is  consistent  with  a  theory  for  scattering  of  radio  waves  by 
field-aligned  ionization  proposed  by  Booker  (ref.  3). 

Two  previous  mathematical  developments  of  forward  scatter  geometry  by  Leadabrand 
and  Yarbroff  (ref.  4)  and  Millman  (refs.  5,  6)  used  the  specular  reflection 
condition.  Both  treatments  were  confined  to  a  single  scattering  point  in  the 
northern  geomagnetic  hemisphere,  Leadabrand  and  Yarbroff  using  plane  earth 
geometry  while  Millman  included  the  curvature  of  the  earth  by  using  a  spherical 
co-ordinate  system. 

In  order  to  apply  the  scatter  geometry  as  a  possible  explanation  of  VHF  radio 
propagation  over  transequatorial  paths,  Millman's  equations  have  been  checked 
for  applicability  to  specular  scattering  in  either  geomagnetic  hemisphere.  In 
addition,  the  present  treatment  includes  paths  which  involve  two  scattering 
irregularities  in  the  ionosphere,  one  on  either  side  of  the  geomagnetic  equator. 
Transequatorial  paths  can  thus  be  calculated  for  radio  waves  which  are  scattered 
at  either  one  or  two  ionospheric  points  before  returning  to  the  surface  of  the 
earth. 
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2.  PREVIOUS  SINGLE  SCATTERING  POINT  GEOMETRY 


The  first  step  is  to  introduce  the  equations  given  by  Millman  (ref.  5)  and 
verify  their  applicability  to  the  southern  geomagnetic  hemisphere.  For  the 
most  part,  the  notation  will  follow  that  of  Millman  or  will  be  a  logical 
extension  of  it. 


Millman  defines  the  angle  between  the  direction  of  the  magnetic  field  line  and 
the  direction  of  electromagnetic  propagation  at  the  scattering  point  as  the 
propagation  angle  9.  This  is  illustrated  in  figure  1  taken  from  Millman  (ref. 
S). 


The  various  symbols  refer  to  the  following: 

Ps  the  scattering  point  in  space 

Pj-  the  point  vertically  beneath  Pg  on  the  earth's  surface 

R  slant  range  from  the  transmitter  location  to  Pg 

I  magnetic  inclination  angle  (or  dip)  at  point  Pg 

a  the  point  of  intersection  of  the  direction  of  electromagnetic 

propagation  with  the  plane  tangent  to  the  earth's  surface  at  PG 

b  the  point  of  intersection  of  the  direction  of  the  magnetic  field  line 
which  passes  through  Pg 

P^-b  line  in  the  tangent  plane  whose  orientation  coincides  with  a  meridian 
of  geomagnetic  longitude 

a  azimuth  bearing  of  the  transmitter  location  with  respect  to 

e  angle  formed  by  the  intersection  of  the  propagation  vector  and  the 
vertical  connecting  Pg  and  PG  (NB:  the  latter  is  a  radial  from  the 
centre  of  the  earth) 

r  radius  of  the  earth, 
o 


Millman  considers  the  two  triangles  aP„b  and  aP„b  to  obtain  from  the  law  of 

b  U 

cosines , 


and 


_ 2  _  2  _  2  _  _ 

ab  =  aPg  +  bPg  -  2  aPg.bPg  cos(180  -  6) 


apG2  ♦  bPG 


2  aPG.bPG  cos  a 


(1) 


ab 


2 


(2) 


-  3  - 
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Equating  the  two  expressions,  it  follows  that 

2  _  _  2  _  _ 

h  +  2  aPg.bPg  cosG  =  -h  -  2  aP^.bP^,  cosa 


where 


h  =  PSPG  aPS  '  aPG  bPS  '  bPG 


Since 


aP, 


sin  e  = 


cos  e  = 


aPr 


aPr 


T  h  bPG 

sin  I  =  -  cos  I  =  - 


bP, 


bP, 


S  S 

equation  (3)  simplifies  to 

cose  =  -  cos  e  sin  I  -  sin  e  cos  I  cosa 


(3) 

(4) 


(5) 


Mi liman  then  specifies  the  angles  e,  I  and  a  in  terms  of  known  parameters  to 
solve  for  e.  The  parameters  he  assumes  are  XR  and  4>R,  the  geographic  longitude 
and  latitude  of  the  transmitter  site  and  the  co-ordinates  of  point  Pg  with 
respect  to  that  location;  A,  E  and  h.  He  derives  the  magnetic  field  quantities 
I  and  a  from  an  earth-centred  dipole  model  of  the  field  with  the  geomagnetic 
north  pole  at  geographic  longitude  A^  and  latitude  i|>M. 

The  inclination,  or  magnetic  dip  I,  specifies  the  direi ' ion  of  the  total  magnetic 
intensity  vector  with  respect  to  the  horizontal  and,  by  definition,  is  measured 
positively  in  the  downward  direction.  For  a  dipole  field  model,  it  is  a  function 
only  of  geomagnetic  latitude  defined  by 

-l 

I  =  tan  (2  tani)>p)  ,  (6) 

where  ipp  is  the  geomagnetic  latitude  of  both  points  Pg  and  PG  and  can  be 
expressed  in  terms  of  geographic  co-ordinates  $p,  Ap. 

-l 

*p  =  sin  {cos(XM  -  Ap)  cos<i>M  cos<j>p  +  sin<f>M  sin$p}  (7) 

In  Millman’s  later  paper  (ref.  6)  he  defines  a  by  the  equation 

a  =  y  -  D  ,  (8) 

where 

y  is  the  geographic  azimuth  (bearing)  of  the  transmitter  location  with 
respect  to  Pfi  and 
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D  is  the  local  magnetic  declination  (or  variation). 

Both  of  these  may  be  expressed  in  terms  of  geographic  co-ordinates 


Y  =  tan 


sin(XR  -  Xp)  cos<j>R 

sin$R  cos<j>p  -  cos<}>R  sin<fp  cos(XR  -  Xp) 


D  =  tan 


sin(XM  -  Xp)  cos<()M 


[sin<f>M  cos^p  -  cos<)>M  sin<t>p  cos(XM  -  Xp) 
the  latter  expression  applying  for  the  dipole  field. 


(9) 


(10) 


Finally,  the  angle  e  can  be  evaluated  by  the  sine  rule  from  the  vertical 
triangle  to  the  earth's  centre. 


e 


-1 

sin 


cos 


E 


(11) 


Having  established  the  incident  propagation  angle,  the  direction  of  propagation 
after  the  forward  scattering  can  be  determined.  As  shown  in  figure  2  taken 
from  reference  6,  Millman  defines  (tt  -  9^)  as  the  propagation  angle  formed  by 
the  scattered  ray  and  the  magnetic-field  line.  Thus  the  condition  of  specular 
scatter  requires  that 


cos0r  =  -  cose  (12) 

3.  APPLICATION  IN  THE  SOUTHERN  HEMISPHERE 

Figure  3  illustrates  the  application  of  the  scatter  geometry  to  a  point  in  the 
southern  geomagnetic  hemisphere.  Some  emphasis  is  placed  on  the  transfer  here, 
since  it  is  not  immediately  obvious  that  the  previous  equations  are  valid  in 
this  situation. 

There  are  two  differences  from  figure  1,  associated  with  the  magnetic  field. 

By  definition,  the  dip  angle  must  be  negative  ir  the  southern  hemisphere.  In 
proceeding  from  equation  (3)  to  equation  (4),  the  absolute  value  of  dip  angle 
was  transferred  into  the  triangle  bPGPg  in  order  to  obtain  expressions  for  the 
functions  sine  and  cosine  of  I.  The  same  transfer  can  again  be  made  but  the 
negative  sign  associated  with  the  dip  angle  must  not  be  forgotten.  The  second 
difference  is  that  the  angle  bP^a  is  no  longer  a,  but  is  equal  to  (it  -  a). 
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A  further  change  is  that  the  angle  at  the  apex  of  triangle  aPgb  becomes  equal  to 
0  rather  than  its  supplement  as  before. 

Thus  the  equation  (3)  is  re-written  as 

2  _  _  2  _  _ 

h  -  2  aPg.bPg  cose  =  -h  +  2  aPG.bPG  coso  ,  (13) 

which  after  substitution  becomes 

cos0  =  cos  e  s in ( I [  -  sin  e  cos|l|cosa  .  (14) 

Since  1  will  be  negative  and  will  lit  in  the  range  -  —  <  I  <  0,  the  equation 
(14)  may  be  converted  to 

cos©  =  -  cos  e  sin  I  -  sin  e  cos  I  cosa  ,  (15) 

which  is  exactly  the  same  as  equation  (5). 


4.  TRANSEQUATOR I AL  SCATTERED  RAY  PATHS 


Two  types  of  transequatorial  paths  are  proposed.  The  first  are  those  families 
of  paths  in  which  the  rays  undergo  forward  scattering  at  a  single  point  and 
return  to  the  earth's  surface  in  the  opposite  geomagnetic  hemisphere  to  the 
transmitter  location.  The  other  type  include  the  rays  which  miss  the  earth 
after  one  forward  scattering,  but  then  reach  ionospheric  heights  in  the  opposite 
geomagnetic  hemisphere.  Under  appropriate  conditions,  forward  scattering  could 
occur  at  a  second  point,  whence  some  of  the  ray  paths  would  return  to  the  earth's 
surface  at  much  higher  latitudes  in  the  second  hemisphere. 


It  is  possible  to  trace  those  ray  paths  which  could  be  transmitted  from  a 
certain  location  <(> _ ,  XD,  and  propagate  via  this  mechanism  to  points  in  the 
opposite  hemisphere  by  considering  ranges  of  initial  elevation  angle  E  and 
geographic  azimuth  bearing  C.  Heights  at  which  the  scattering  might  occur  are 
fixed  at  h^  in  the  transmitting  hemisphere  and  h ^  in  the  other  hemisphere.  The 
angle  e  between  the  initial  ray  path  and  the  zenith  at  the  first  point  of 
forward  scattering  is  determined  from  equation  (11),  i.e. 


e 


cos 


E 


(16) 
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It  follows  that  the  earth  centre  angle  between  the  transmitter  location  and  the 
scattering  point  is 

v  =  j  -  E  -  e  .  (17) 

The  co-ordinates  of  the  first  scattering  point  can  then  be  obtained  by  applying 
the  sine  and  cosine  rules  for  spherical  triangles,  viz 


and 


-1 


=  sin 


sim)>  cosv  +  cos4>  sinv  cos£ 
R  R 


X 


p 


X^  +  sin 


sin£  sinv 
cos$p 


(18) 


(19) 


The  bearing  of  the  transmitter  location  with  respect  to  the  sub- ionospheric 
point,  i.e.  the  location  on  the  earth's  surface  directly  beneath  the  first 
scattering  point,  will  be  y  as  given  by  equation  (9). 


Thus  the  incident  propagation  angle  0  can  be  found  using  equation  (3),  allowing 
substitution  of  equation  (8). 


0 


-l 

cos 


' 

-cos  e  sin  1  -  sin  e  cos  I  cos(y 

. 


D) 


(20) 


It  remains  to  specify  the  magnetic  field  parameters  I  and  D,  which  may  be  done 
simply  by  the  dipole  field  model  as  given  before,  using  particular  geographic 
co-ordinates  for  the  north  geomagnetic  pole,  thus 

Xw  =  70.1°  xl 
M 

<t>M  =  78.6°  M  .  (21) 

Alternatively,  the  magnetic  field  parameter;  could  be  derived  during  any 
computation  from  a  spherical  harmonic  repre; entation  of  the  field  such  as  that 
of  Jensen  and  Cain  (ref.  7).  It  would  also  be  possible  to  determine  these 
parameters  by  interpolation  of  values  scalec  from  magnetic  charts  over  a 
suita >le  grid,  having  stored  the  grid  values  in  a  matrix  for  computer  reference. 

Once  the  forward  scattering  has  occurred,  the  scattered  signal  described  a 
conical  surface  surrounding  the  magnetic  field  line,  with  the  apex  at  the 
scattering  point.  All  the  scattered  rays  in  the  cone  have  the  same  propagation 
angle  of  (it  -  6  )  as  defined  in  figure  2.  Depending  on  the  magnetic  dip  angle. 
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some  of  the  lower  part  of  the  cone  might  intersect  the  surface  of  the  earth  in  a 
smooth  curve,  while  other  ray  paths  can  miss  the  earth  and  propagate  into  the 
opposite  hemisphere. 

The  expression  for  the  scattered  propagation  angle  is  similar  to  equation  (20) 

cos(u  -  8^)  =  -  cos  er  sin  I  -  sin  er  cos  I  cos(vr  -  D)  ,  (22) 

where  the  subscript  r  refers  to  the  scattered  (or  reflected)  ray. 

The  equation  may  be  re-arranged  to  obtain  y  ,  which  is  the  azimuth  of  the 
scattered  rays 


Y  =  cos 
r 


cos6  -  cos  e  sin  I 
_ r _ r _ 

sin  e_  cos  I 
r 


Substituting  equation  (12), 


-i  -  cos8  -  cos  e_  sin  I 

1  T 


Y  =  cos 
r 


sin  er  cos  I 


The  cone  is  formed  by  allowing  the  zenith  angle  er  to  vary  over  a  range  of 
values  from 


<*r>.in  =  2  '  (I  +  9) 


(e  )  =  (e  )  .  +  28  .  (26) 

v  r-'max  r;min 

Both  these  extremes  apply  to  scattered  rays  lying  in  the  direction  of  the 
magnetic  field,  vertically  below  and  above  the  field  line.  The  inverse  cosine 
function  of  equation  (24)  may  take  either  a  positive  or  negative  value  to 
produce  the  horizontal  spread  of  the  cone  on  each  side  of  the  field  line. 

The  scattered  rays  which  intersect  the  earth's  surface  must  lie  below 


(e  ).  =  sin  - r— 

r  t  v0  +  hj 


which  is  the  condition  of  tangency  with  the  earth. 


If  (e  )  .  is  less  than  (e  )„,  then  part  of  the  cone  will  intersect  the  earth's 
v  ruin  v  r-’t’  r 

surface  and  the  location  at  which  the  rays  reach  the  surface  may  be  calculated. 
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The  earth  centre  angle  between  the  scattering  point  and  each  earth  intersection 
point  is  given  by  g  . 


6r  =  sin 


r  +  h, 
o  1 

-  sin  e 

r  r 


-  e 


(28) 


The  geographic  co-ordinates  of  these  points  are  then 

'll 


=  sm 


sin^p  cos$r  +  cos<j>p  sing^  cosy^ 


.  (29) 


-1 


X  =  Xn  +  sin 
r  P 


sing  siny 
r  r 


COS<t> 


(30) 


The  scattered  rays  in  the  cone,  which  will  be  forward  scattered  a  second  time 
at  the  appropriate  height  in  the  opposite  hemisphere,  can  now  be  considered. 
Figure  4  introduces  the  notation  to  be  used  by  illustrating  a  transequatorial 
ray  path  lying  in  the  geomagnetic  longitude  meridian  plane  of  the  simple  dipole 
field  model.  Similar  subscripts  are  used  but  the  prime  (')  indicates  relation 
to  the  second  scattering  point. 


The  second  scatter  point  is  determined  by  the  triangle  with  its  base  from  the 
first  to  the  second  scatter  point  and  its  apex  at  the  centre  of  the  earth.  It 
is  thus  possible  to  transpose  to  the  second  scatter  point  using 


e ' 


sin 


r  +  h, 
o  1 

- r —  sin  e 

r  +  h-  r 


(31) 


From  the  same  triangle,  the  earth  centre  angle  between  the  two  scatter  points  is 


er  '  e 


(32) 


and  bearing  from  the  first  to  the  second  is  yr  as  given  in  equation  (24).  Thus 

'K 


=  sin 


s:n<}>p  cosu)  +  cos4>p  sinu  cosyr 


(33) 


and 


»p  =  A|>  +  sin 


-1 


sinw  siny 


cos4>  ’ 
P 


(34) 


The  range  of  values  of  e^  which  will  produce  such  ray  paths  begins  at  the 
tangency  condition  of  equation  (27),  with  the  upper  limit  occurring  when  the 
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ray  path  reaches  the  second  scattering  height  at  the  geomagnetic  equator. 

If  the  dipole  field  relation 

~l 

ijjp  =  tan  (h  tan  I)  (35) 

is  assumed  for  any  field  model,  the  upper  limit  can  be  expressed  by  the  implicit 
relation 


ro  ♦  hi 

~nr  Sin  er  sin(V  +  er)  •  (36) 

o  2 


A  number  of  second  scatter  points  will  be  found  and,  for  each  location  C <t>p ,  Xp) , 
the  appropriate  magnetic  field  parameters  determined.  These  are  substituted  in 
an  expression  for  the  second  incident  propagation  angle  0'. 

-U  ) 


9 '  =  cos 


-  cos  e'  sin  I'  -  sin  e'  cos  I'  cos(y'  -  D') 


(37) 


where  D'  is  the  magnetic  declination  at  the  second  scatter  point  and 

Y '  is  the  bearing  of  the  first  scatter  point  with  respect  to  the  second. 


Y ' 


-1 

tan 


sin(X  -  X')  cost|> 

_ .11.  .£  ...2 _ 

s in<t>p  cosi^p  -  cos^p  sin^  cos(Xp  -  X£) 


(38) 


The  azimuth  bearing  of  the  ray  paths  after  a  second  forward  scattering  is  given 
by 


-1 


cos 


-  cos8  ' 


cos  e '  sin  I ' 
r 


sin  e  ’  cos  I ' 
r 


+  D' 


(39) 


Therefore,  at  each  of  the  possible  second  scatter  points  new  cones  of  rays  will 
be  produced,  some  of  which  will  intersect  the  earth's  surface  to  form  sets  of 
smooth  arcs. 


The  limits  on  values  e^  may  take  are 


(e')  . 
v  rJmin 


-  -  (6'  -  I') 


and 


<e;>t 


-i 

sin 


(40) 


(41) 
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Finally,  the  sets  of  earth  centre  angles  between  the  second  scatter  points  and 
the  earth  intersection  points  of  their  conical  scatter  surfaces  are 


8'  = 
r 


sm 


♦  h. 


sin  e  1 


e  * 
r 


(42) 


which  leads  to  the  general  expressions  for  the  geographic  co-ordinates  of  these 
intersect  ions 


and 


t> '  =  sin 
r 


sin4>^  cos8^  +  cos^  sinB^  cosy^ 


X ' 
r 


-l 

+  sin 


sing'  siny' 
r  r 

cos<f> ' 
r 


(43) 


(44) 


5.  CONCLUDING  REMARKS 

The  mathematical  description  of  field-aligned  ionization  scatter  geometry, 
originally  derived  by  Millman  (refs.  5,  6)  for  forward  scattering  of  radio  waves, 
has  been  extended  to  include  transequatorial  paths. 

Mathematical  expressions  have  been  derived  which  will  enable  the  calculation  of 
ray  paths  from  a  transmitter  in  one  geomagnetic  hemisphere  to  points  on  the 
surface  of  the  earth  in  the  other  hemisphere,  where  the  mode  of  propagation  has 
been  by  means  of  forward  specular  scattering. 

There  can  be  two  types  of  transequatorial  scattered  ray  paths.  The  first  undergo 
forward  scattering  at  a  point  in  the  same  hemisphere  as  the  transmitter,  but 
return  to  ground  as  the  lower  portion  of  the  scattered  ray  cone,  the  geographic 
co-ordinates  of  the  resulting  arc  of  intersection  being  given  by  equations  (29) 
and  (30) .  The  other  type  of  path  arises  from  the  two  arcs  of  the  same  scatter 
cone  horizontally  on  each  side  of  the  central  field  line,  which  miss  the  earth, 
but  intersect  the  ionosphere  in  the  opposite  hemisphere  to  produce  a  set  of 
second  scattering  points.  Each  of  these  in  turn  results  in  a  cone  of  scattered 
ray  paths,  part  of  which  may  intersect  the  earth's  surface.  The  locii  of  the 
intersection  points  after  the  secondary  forward  scatter  are  sets  of  smooth  arcs, 
the  geographic  co-ordinates  of  which  can  be  obtained  from  equations  (43)  and  (44). 

It  should  be  possible  to  trace  all  such  ray  paths  numerically,  employing  a 
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digital  computer.  Subsequent  analysis  would  establish  whether  the  mechanism 
offers  a  feasible  explanation  of  certain  anomalous  VHF  transequatorial 
propagation  which  has  been  observed  over  experimental  paths. 
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